Nonmonotonic temperature dependence of critical current in diffusive d-wave 

junctions 
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We study the Josephson effect in D/I/DN/I/D junctions, where I, DN and D denote an insu- 
lator, a diffusive normal metal and a d-wave superconductor, respectively. The Josephson current 
is calculated based on the quasiclassical Green's function theory with a general boundary condi- 
tion for unconventional superconducting junctions. In contrast to s-wave junctions, the product 
of the Josephson current and the normal state resistance is enhanced by making the interface bar- 
riers stronger. The Josephson current has a nonmonotonic temperature dependence due to the 
competition between the proximity effect and the midgap Andreev resonant states. 



Since the discovery of Josephson effect 1 in supercon- 
ductor / insulator / superconductor (SIS) junctions, it 
has been studied in various types of structures^. The 
Josephson current flows not only through thin insulators 
but also through normal metals or ferromagnets^* 3 ^. In 
superconductor / diffusive normal metal / superconduc- 
tor (S/DN/S) junctions, Josephson current is carried by 
Cooper pairs penetrating into the DN as a result of the 
proximity effect. In s-wave superconductor junctions, the 
maximum amplitude of Josephson current (Ic) mono- 
tonically increases with the decrease of temperatur o 5 i 6 i 7 . 
This is a natural consequence of the fact that interference 
effects arc stronger at lower temperatures. The Joseph- 
son effect, however, depends strongly on pairing sym- 
metries of superconductors because it is an essentially 
phase sensitive phenomenon. In d-wave superconductor 
/ insulator / d-wave superconductor (DID) junctions, for 
instance, Ic first increases with the decrease of temper- 
ature then decreases and can even change its sign at 
low temperatures for certain oricntations&^iiLII. This 
nonmonotonic behavior of Josephson current is caused 
by 0-7T transition due to the mid-gap Andreev resonant 
states (MARS) forming at junction interfaces'^ 3 -. Sim- 
ilar effect is also observed in Josephson current through 
ferromagnets*. Thus the nonmonotonic temperature de- 
pendence of Ic is a sign of unusual interference effect. 

The quasiclassical Green's function theory is a power- 
ful tool to study quantum transport in superconducting 
junctions. The circuit theory^ provides the boundary 
conditions for the Usadel equations^ widely used in dif- 
fusive superconducting junctions. These boundary condi- 
tions generalize the Kupriyanov-Lukichev conditions 6 for 
an arbitrary type of connector which couples the diffusive 
nodes. Recently Tanaka et almiiJ^i— have extended the 
circuit theory^ to systems with unconventional super- 
conductors. An application of the extended circuit the- 
ory to DN/d-wave superconductor junctions has revealed 
a strong competition (destructive interference) between 
the MARS and the proximity effect in DNlSJi. This 
competition, however, has not yet tested experimentally. 
Thus it's important to propose the way to verify this 



prediction. 

In the present paper, we show that Josephson effect is 
a suitable tool to observe the above competition. We ex- 
tend the previous theoryA&i£ in order to treat the exter- 
nal phase difference between the left and right supercon- 
ductors. Applying this formalism, we calculate Joseph- 
son current in d-wave superconductor /insulator/ dif- 
fusive normal metal /insulator/ d-wave superconductor 
(D/I/DN/I/D) junctions, solving the Usadel equations 
with new boundary conditions. This allows us to study 
the influences of the proximity effect and the formation 
of MARS on the Josephson current simultaneously. We 
find that the competition between the proximity effect 
and the formation of MARS provides a new mechanism 
for a nonmonotonic temperature dependence of the max- 
imum Josephson current. 

We consider a junction consisting of d-wave super- 
conducting reservoirs (D) connected by a quasi-one- 
dimensional diffusive conductor (DN) with a resistance 
Rd and a length L much larger than the mean free path. 
The DN/D interface located at x = has the resistance 
R' b , while the DN/D interface located at x = L has the 
resistance Rb- We model infinitely narrow insulating bar- 
riers by the delta function U(x) = HS(x — L) + H'8(x). 
The resulting transparencies of the junctions T m and 
T' m are given by T m = 4 cos 2 </>/(4cos 2 + Z 2 ) and 
T' m = 4cos 2 (/)/(4cos 2 + Z' 2 ), where Z = 2H/v F and 
Z' = 2H'/vjr are dimensionless constants and <f> is the in- 
jection angle measured from the interface normal to the 
junction and vf is Fermi velocity. We show the schematic 
illustration of the model in Fig. 2] 

We parameterize the quasiclassical Green's functions 
G and F with a function <&J&: 
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where u) is the Matsubara frequency. Then Usadel equa- 
tion reads^ 
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FIG. 1: (color online) Schematic illustration of the model. 



with the coherence length £ = ^JD/V/kTc, the diffusion 
constant D and the transition temperature Tq- From the 
conservation law for the matrix currenljii, the boundary 
conditions are given by 
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in a time reversal symmetry conserving state. Here, S z 
denotes the z-component of the total spin of a Cooper 
pair. 

The average over the various angles of injected parti- 
cles at the interface is defined as 
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with B(<j)) = B and T (,) (cj)) = T$. It is important to 
note that the solution of the Usadel equation is invariant 
under the transformation a — > —a or — * —j3. This is 
clear by replacing <f> with —(f) in the angular averaging. 
The resistance of the interface is given by 
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denotes R b or R' b , 



and R^ is 



Here, for example, 
Sharvin resistance, which in three-dimensional case is 
given by R^~ l = e 2 k 2 F Sc^ /(47r 2 ), where kp is the Fermi 
wave-vector and Sc is the constriction area. 
The Josephson current is given by 
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g ± = w/y/tf + Al f± = A±/ 
A(T) cos(2^ =F 2a) at x = and 
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at x = L, where Ii,l2 and 7 3 are defined similarly to 
I[, I 2 and ^3 by removing the superscript "', exchanging 
subscript '+' for subscript'-', putting tp — and substi- 
tuting f3 for a, respectively. Here ip is the external phase 
difference across the junctions, and a and f3 denote the 
angles between the normal to the interface and the crys- 
tal axes of eZ-wave superconductors for x < and x > L 
respectively. 

These boundary conditions derived above are quite 
general since with a proper choice of A± they are appli- 
cable to any unconventional superconductor with S z = 
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with temperature T and R = Rd + Rb + R' b - In the 
following we focus on the IcR value where Ic denotes 
the magnitude of the maximum Josephson current. Wc 
consider symmetric junctions with R b = R' b and Z = Z' 
for simplicity. 

In Fig. HJa), IcR is plotted as a function of tempera- 
ture for R d /R b = 1, E Th /A(0) = 1, and {a, S3) = (0,0). 
The results show that IcR increases with the decrease of 
T as it does in s-wave junctions as shown in (b). Am- 
plitude of IcR increases with the increase of Z in d- 
wave junctions, whereas the opposite tendency can be 
seen in s-wave junctions. The sign change of pair poten- 
tial is responsible for the Z dependence of IcR in (a). 
In the d-wave symmetry with a = = 0, injection an- 
gles of a quasiparticle can be divided into two regions: 
4>+ = {>|0 < |0| < tt/4} and 0_ = {0|tt/4 < \<f>\ < tt/2}. 
The sign of pair potential for (f> + and that for are 
opposite to each other. In general, such a sign change of 
pair potentials leads to the suppression of the proximity 
effect in DN and hence Josephson currents. This is be- 
cause Cooper pairs originated from the positive part of 
pair potentials cancel those originated from the negative 
part of pair potentials. For small Z, a quasiparticle with 
0+ and that with (/)_ give opposite contribution to the 
proximity effect. In the actual calculation, the cancella- 
tion due to the sign change of the pair potential reduces 
the magnitudes of I 2 , I3, I 2 , and I' 3 . On the other hand 
for large Z, integration over has the dominant contri- 
bution to I 2 , I3, I 2 , and I' 3 . As a result, the cancellation 
of Cooper pairs becomes negligible in low transparent 
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junctions. Thus IcR in the d-wave junctions increases 
with increasing Z . This argument can be confirmed by 
the behavior of $ w which represents a degree of the prox- 
imity effect in DN. In Figs.|2c) and (d), we show $ w as 
a function of cj, where we choose x — L/2, ip = tt/2, 
T/T c = 0.2, and (a, (3) = (0,0). The magnitudes of 
Re $ w and im<I> w increase with increasing Z. This is 
the reason of the enhancement of IcR product for low 
transparent junctions with large magnitude of Z . 
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FIG. 2: (color online) (a) Temperature dependence of IcR 
for (a, (3) = (0,0). (b) The results for s-wave junctions for 
comparison. Real (c) and imaginary (d) parts of $^ as a 
function of u>. Here we choose Rd/Rb = 1 and Eth/A.(0) = 1. 

Next we focus on the Josephson effect with nonzero val- 
ues of a and (3. Figure|21(a) displays IcR as a function of 
temperature for Z = 10, Rd/Rb — 1 and Erfc/A(0) = 1- 
In contrast to IcR for (a,(3) — (0,0), the results for 
(a,/?) = (7r/8, 0) and (n/8,n/8) show a nonmonotonic 
temperature dependence. The transparency at the in- 
terfaces greatly changes the peak structure as shown in 
(b), where temperature dependence of IcR is plotted for 
several Z at (a,/3) = (tt/8,0). The peak structure is 
smeared as the decrease of Z. 

To understand the origin of the nonmonotonic temper- 
ature dependence of IcR, we have to consider in detail 
the mechanism of competition between the proximity ef- 
fect and the formation of the MARS. For example, at a = 
(3 = 7r/8, the injection angles of a quasiparticle are sep- 
arated into four regions: (i) <pM — {4*^/8 < 4> < 37r/8}, 
(ii) cf)p = {(j)\0 < cf) < tt/8 and 3tt/8 < <j> < tt/2}, (iii) 
—(/>m, and (iv) —<f>p. Transport channels within ±4>p 
contribute to the proximity effect but do not form the 
MARS. On the other hand, transport channels within 
±4>m are responsible for the formation of the MARS but 
do not contribute to the proximity effect. Thus channels 
within ±(f>p and those within ±0jvf play different roles 
in electron transport. At the same time, the types of 
wave transmission are different for these two transport 
channels. The channels within ±0p are open for the 
normal transmission. On the other hand, in the chan- 
nels within ±4>m, the resonant transmission is also pos- 
sible for small u>. At high temperatures, effects of the 
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FIG. 3: (color online) Temperature dependence of (a) IcR at 
Z — 10 for several orientation angles, (b) IcR at (a, (3) = 
(71/8, 0) for several Z values. Real and imaginary parts of 
as a function of u for (a, (3) — (tt/8,0) and ip = tt/2 at 
x = L/2 are shown in (c) and (d) respectively. Here we choose 
R d /R b = 1 and E Th /A(0) = 1. 



MARS are negligible and Josephson currents are carried 
through channels within ±</>p. At low temperatures, the 
resonant transmission through ±4>m governs electric cur- 
rents. This is a result of an important property of res- 
onant transport. When a normal transport channel and 
a resonant transport channel are available in parallel, a 
quasiparticle tends to choose the resonant channel for its 
transmission at low temperatures. This property can be 
confirmed in the present junctions by u> dependence of 
$ as shown in Fig.[3](c) and (d). The results show that 
<& rapidly decreases to zero in the limit of u) = 0, which 
implies no diffusion of Cooper pairs into DN through the 
normal channels within ±rfip. At u> = 0, quasiparticles 
choose the resonant channels within ±0m- Electric cur- 
rents through i(j)M, however, do not contribute to the 
net Josephson current because of cancellation of the cur- 
rents through ipM by those through —</>m- Thus $ goes 
to zero at lo = 0, which leads to the suppression of the 
Josephson current at low temperatures and hence the 
nonmonotonic temperature dependence of Ic- The in- 
fluence of the MARS is more remarkable for larger Z. 
Therefore the resulting nonmonotonic temperature de- 
pendence becomes much more pronounced for large Z as 
shown in Fig. CUb) . The origin of the nonmonotonic tem- 
perature dependence of IcR in D/I/DN/I/D junctions 
is different from that in the mirror-type DID junctions 
(i.e., a = — j3), where the current through the resonant 
channels causes the 0-7T transition and the nonmonotonic 
temperature dependenca&IfiiAi. 

It is important to note that the angle average in Eq. (J5J 
at the two interfaces can be carried out independently 
The physical meaning can be explained by considering 
the motion of a quasiparticle which starts from the left 
interface with a certain injection angle and reaches the 
right interface after traveling across the DN. At the right 
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interface, injection angles of a quasiparticle are indepen- 
dent of the initial injection angle at the left interface be- 
cause of the diffusive scatterings in the DN. This property 
allows one to carry out the angular averaging at the two 
interfaces independently. When a DN is replaced by a 
clean normal metal or a clean insulator, the injection an- 
gles at both interfaces are correlated to each other. This 
is because the injection angle of a quasiparticle at the left 
interface is conserved when the quasiparticle reaches the 
right interface. As a result, in the clean limit, Josephson 
current in a symmetric junction (i.e., a — (3) has a be- 
havior different from that in a mirror-type junction with 
a = —(3. In a symmetric DID junction, IcR increases 
monotonically with decreasing temperature^. In the 
D/I/DN/I/D junctions considered in the present paper, 
the mirror-type and the symmetric configurations yield 
the same Josephson current because of the independent 
angle averaging at the two interfaces. 

To study the nonmonotonic temperature dependence 
of IcR further, we define the temperature T p at which 
IcR has a maximum. We have confirmed that T p is in- 
sensitive to Z when Z is sufficiently large. The calcu- 
lated T p shown in Fig. 01 increases monotonically with 
the increase of Rd/Rb and EthI 'A(Q). These results im- 
ply that large magnitudes of Rd/Rb and £'t/i/A(0) 
are needed to observe the peak structure experimentally. 
These conditions are satisfied in junctions with high in- 
sulating barriers, small cross section area of the DN and 
thin DN. 
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FIG. 4: The peak temperature as a function of Rd/Rb (a) and 
£W A (0) (b) for Z = 10 and (a,/3) = (tt/8,0). 



In the present paper, we have studied the Josephson ef- 
fect in d-wave superconductor / insulator / diffusive nor- 
mal metal / insulator / d -wave superconductor junctions 
by solving the Usadel equations. We have derived the 
boundary conditions for the quasiclassical Green's func- 
tion which make it possible to calculate the Josephson 
current across the structure. By calculating the Joseph- 
son current for various parameters, we have clarified the 
following points. (1) In contrast to the conventional s- 
wave junctions, the magnitude of IcR is enhanced with 
the increase of the insulating barrier strengths at the in- 
terfaces. (2) The IcR value may exhibit a nonmono- 
tonic temperature dependence due to the competition 
between the proximity effect and the formation of the 
MARS. The origin of such a nonmonotonic behavior is 
completely different from that in the established case of 
a clean DID junctior&i2*ii where 0-ir transition occurs 
due to the MARS. In order to observe experimentally the 
predicted nonmonotonic temperature dependence of Ic 
in D/I/DN/I/D junctions, large magnitudes of parame- 
ters Z, Rd/Rb, and -&r/i/A(0) are required. Note that we 
can reproduce the main conclusions obtained here within 
the quasiclassical theory of superconductivity by numer- 
ical simulations based on the recursive Green function 
methodic. 

It follows from the present approach that similar non- 
monotonic temperature dependence is also expected for 
s-wave superconductor / insulator / diffusive normal 
metal / insulator / d -wave superconductor junctions. To 
the best of our knowledge of Josephson junctions consist- 
ing of s- and d -wave superconductors, there is only one 
report in literature about the observation of nonmono- 
tonic temperature dependence of Ic in YBaCuO/Pb 
junctions^. Though this experiment might be possi- 
bly relevant to our results, more experimental data are 
needed to check whether it is the case. 



1 B. D. Josephson, Phys. Lett. 1, 251 (1962). 

2 K.K. Likharev, Rev. Mod. Phys. 51, 101 (1979). 

3 A. A. Golubov, M. Yu. Kupriyanov, and E. Il'ichev Rev. 
Mod. Phys. 76, 411 (2004). 

4 V. V. Ryazanov, V. A. Oboznov, A. Yu. Rusanov, A. V. 
Veretennikov, A. A. Golubov, and J. Aarts, Phys Rev. 
Lett. 86, 2427 (2001). 

5 V. Ambegaokar and A. Baratoff, Phys. Rev. Lett. 10, 486 
(1963). 

6 M. Yu. Kupriyanov and V. F. Lukichev, Zh. Exp. Teor. 
Fiz. 94 (1988) 139 [Sov. Phys. JETP 67, (1988) 1163]. 

7 A. V. Zaitsev, Physica C 185-189, 2539 (1991). 

8 Yu. S. Barash, H. Burkhardt, and D. Rainer, Phys. Rev. 
Lett. 77, 4070 (1996); Y. Tanaka and S. Kashiwaya, Phys. 
Rev. B 53, R11957 (1996); Y. Tanaka and S. Kashiwaya, 
Phys. Rev. B 56, 892 (1997). 



9 A. A. Golubov and M. Yu. Kupriyanov, JETP Lett. 69, 
262 (1999). 

E. Il'ichev, M. Grajcar, R. Hlubina, et ai, Phys. Rev. Lett. 
86, 5369 (2001). 

1 G. Testa, E. Sarnelli, A. Monaco, E. Esposito, M. Ejr- 
naes, D.-J. Kang, S. H. Mennema, E. J. Tarte, and M. G. 
Blamire Phys. Rev. B 71, 134520 (2005) 

2 L.J. Buchholtz and G. Zwicknagl, Phys. Rev. B 23 5788 
(1981); J. Hara and K. Nagai, Prog. Theor. Phys. 74 
(1986) 1237 ; C. Bruder, Phys. Rev. B 41, (1990) 4017; 
CR. Hu, Phys. Rev. Lett. 72, (1994) 1526. 

3 Y. Tanaka and S. Kashiwaya, Phys. Rev. Lett. 74, (1995) 
3451. 

4 Yu. V. Nazarov, Superlattices and Microstructuctures 25, 
1221 (1999). 

5 K.D. Usadel, Phys. Rev. Lett. 25, 507 (1970). 



■5 



Y. Tanaka, Y.V. Nazarov and S. Kashiwaya, Phys. Rev. 
Lett. 90, 167003 (2003). 

Y. Tanaka, Yu. V. Nazarov, A. A. Golubov, and S. Kashi- 
waya, Phys. Rev. B 69, 144519 (2004). 
Y. Tanaka and S. Kashiwaya, Phys. Rev. B 70, 012507 
(2004); Y. Tanaka, S. Kashiwaya and T. Yokoyama, Phys. 
Rev. B 71, 094513 (2005); Y. Tanaka, Y. Asano, A. A. 



Golubov and S. Kashiwaya, Phys. Rev. B 72, 140503(R) 
(2005). 

Y. Asano, Phys. Rev. B 64, 014511 (2001); J. Phys. Soc. 
Jpn 71, 905 (2002). 

I. Iguchi and Z. Wen, Phys. Rev. B 49, R12388 (1994). 



